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We investigate the critical transition from an inverse cascade of energy to a forward energy 
cascade in a two-dimensional magneto-hydrodynamic flow as the ratio of magnetic to mechanical 
forcing amplitude is varied. It is found that the critical transition is the result of two competing 
processes. The first process is due to hydrodynamic interactions, cascades the energy to the large 
scales. The second process couples small scale magnetic fields to large scale flows transferring the 
energy back to the small scales via a non-local mechanism. At marginality the two cascades are 
both present and cancel each other. The phase space diagram of the transition is sketched. 

PACS numbers: May be entered using the \pacs{#lj command. 


I. INTRODUCTION 

The cascade of ideal invariants across scales, is a fun¬ 
damental concept of turbulence theory. In three dimen¬ 
sional flows energy and helicity cascade forward to the 
small scales while in two dimensions energy cascades in¬ 
versely to the large scales and enstrophy cascades for¬ 
ward to the small scales. Many flows in nature how¬ 
ever [1] show both characteristics: the formation of large 
structures due to an inverse cascade of energy and small 
scale turbulence due to a forward energy cascade. Coex¬ 
istence of forward and inverse cascades has been observed 
in laboratory and in numerical experiments under differ¬ 
ent physical situations: confined turbulence in thin lay¬ 
ers, rotating and stratified turbulence, turbulence in the 
presence of strong magnetic fields, 2D turbulence in the 
presence of magnetic forcing [m. In these studies the 
amplitude of the forward and inverse cascade can be var¬ 
ied by changing the relevant control parameter (rotation, 
stratification, geometric factor, magnetic field strength 
etc). Due to high computational costs of three dimen¬ 
sional numerical simulations, it is difficult to determine 
precisely the way the cascade transitions from forward to 
inverse. A notable exception is the case of transition from 
two-dimensional hydrodynamic turbulence (2D HD) to 
two-dimensional magneto-hydrodynamic turbulence (2D 
MHD) [ 3 . In this case the transition is not “dimensional” 
and all the simulations can be carried out in two dimen¬ 
sions which is computationally less expensive. For weak 
magnetic forcing the flow behaves as 2D HD cascading 
energy inversely. As the magnetic forcing is increased 
the inverse cascade decreases and finally stops at a crit¬ 
ical amplitude of the control parameter. Close to this 
critical point the flux of the inverse cascade scales as a 
power law with the distance from criticality. This prop¬ 
erty manifests itself in the limit of large box size while at 
moderate box-sizes the transition appears smooth. 

In this work we unravel the mechanisms involved in 
this transition by looking in detail the different processes 
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involved in transferring ideal invariants across scales both 
in spectral and real space. In the next section 0 
we describe in detail our system and define the non- 
dimensional control parameters and observables that we 
are using in our analysis. Section (HII presents the be¬ 
havior of global (space and time aver aged) quantities 
close to the critical point. Section (IV) shows how the 
fluxes and spectra vary as the control parameter is var¬ 
ied. Section (|^ discusses the effect of the transition on 
spatial structures and their statistics. Finally, in the last 
section we summarize this work, sketch the phase space 
diagram and draw our conclusions. 


II. DESCRIPTION OF THE SYSTEM 

The 2D MHD system is governed by the Navier-Stokes 
equation for the velocity field u and the induction equa¬ 
tion for the magnetic field b. In terms of the vorticity 
w = • V X u and the vector potential a, b = V (xae^), 

the governing equations read: 

dtuj + u- Vw =b • -I- 

dta + u-Wa= -b?7+V2"a-b r?" -b (1) 

where j = • V x b is the current, , 0 ^ are small-scale 

dissipation coefficients while v ~, rj~ are large-scale dissi¬ 
pation coefficients. The parameter n, m give the order of 
the laplacian used in the dissipation terms. The physi¬ 
cally motivated values are n = 1 and m = 0. However, 
in the present work we use hyperviscosity to increase the 
inertial range and fix these values at n = m = 2. These 
equations are numerically solved on a doubly periodic 
2'kL X 2ttL domain using a pseudospectral method. A 
standard Runge-Kutta of fourth order scheme is used for 
time marching (see m for further details on the code). 


A. Nondimensional parameters 

The main control parameter in the present study is 
the ratio of the forcing in the magnetic field to the 
forcing in the velocity field fj,f = Fu is the 
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mechanical forcing given by r„ = —A ^ V x {(/fuiGz) 
and Fb is the magnetic forcing given by Fb = V x 
with 4>^ = 2fQkfCOs{kfx)cos{kfy) and 4>a = 
^ffokJ^sm{kfx)sm{kfy). kf L gives the wave-number 
at which the system is forced, non-dimensionalized by 
the box size. The ratio of the dissipation coefficients re¬ 
sults in two Prandtl numbers: one for the large scale 
Pm~ = 'r]~ I v~ and the one for small scales Pm~^ = 
T]^jv^. Both of these Prandtl numbers are set to 
unity PmP = Pm~ = 1. The strength of turbu¬ 
lence is measured by the Reynolds numbers at small 
and large scales i?e+ = Re“ = 

l\v~\ defined here based on the forcing 

amplitude. ReP' determines the extend of the forward 
energy cascade while Re~ determines the extend of the 
inverse energy cascade. In all runs fc/L was chosen suf¬ 
ficiently large, so that the cascade does not reach the 
box-size and no large scale condensate is formed [IIHS]. 

An alternative control parameter to fif is the ratio of 
the injection rates dehned as ratio of injection energy 
in u to the injection energy in b, /Tg = R/Iu with R = 
(Fb • b) and R = (Fu • u) where the angular brackets 
stand for spatial and time average. This parameter might 
be more fitting to compare with more theoretical models 
like shell/EDQNM models or for forcing functions for 
which the energy injection rates is fixed. However for the 
forcing chosen in this work is not a control parameter 
in our system but an observable. 


B. Observables 

In the ideal 2D hydrodynamic flow there are two 
conserved quantities in the system the kinetic energy 
Eu = i (|up)^ and the enstrophy H ^ (where 

(•)^ stands for spatial average). Eu cascades to larger 
scales and H cascades to smaller scales M- On the other 
hand in the ideal 2D MHD system the two conserved 
quantities are the total energy E = ^ (|up -I- |bp) that 
cascades to the small scales and the square vector poten¬ 
tial A = i that cascades to the large scales [nmi]. 

When forcing and diffusion is included in the system in 
the long time limit the system reaches a steady state 
where the injection rate of all ideally-conserved quanti¬ 
ties is balanced by their dissipation rate in the small and 
large scales. In this case choosing y,/ = 0 reduces the sys¬ 
tem in the long time limit to a 2D HD flow since any ini¬ 
tial magnetic field will disappear due to the anti-dynamo 
theorem for two-dimensional flows m- Accordingly en¬ 
ergy will be dissipated in the large scales while enstro¬ 
phy will be dissipated in the small scales. For non-zero 
values ot yj the magnetic held will be sustained. How¬ 
ever, if the magnetic forcing is sufficiently small yj 1 
the magnetic held will be too weak to feed back on the 
how. In this limit the vector potential acts like a passive 
scalar advected by the 2D HD how. In addition to Ejj 


and H the advection term also conserves the square vec¬ 
tor potential A that cascades to the small scales. As the 
control parameter yf becomes larger, the Lorentz force 
eventually acts back to the how and magnetic held stops 
being passive. Nonlinearities no longer conserve H and 
Ejj since the Lorentz force can inject/absorb kinetic en¬ 
ergy and enstrophy to/from the how. The nonlinearities 
conserve the total energy E and the square vector poten¬ 
tial A. At sufficiently large yf the system transitions to 
2D MHD state that cascades the total energy E forward 
to the small scales while A cascades inversely to the large 
scales. 

The strength of an inverse or a forward cascade of a 
quantity at steady state can be measured by the rate of 
dissipation in the large and small scales respectively. The 
rate of energy dissipation at large and small scales de¬ 
noted by e“, respectively are dehned in Fourier space 
as, 

4^ k+l(Ek^olkp"(|iikp + |bkp))^ 

k-|(Ek^olk|-^™(|akP + |bkp|))^. (2) 

where (•)^ stands for time average. The rate of dissi¬ 
pation of the square vector potential in small and large 
scales denoted by e+, e~ are dehned as 

4^ l^+l(i:k^o|kP”|akp)^ 

£:^k-|(Ek^o|k|-^"‘|akp)^ (3) 

and similarly for the enstrophy dissipation rates we 
dehne 

4^ l^+l(Ek^olkp"+^|Oikp)^ 

e”= I (Ek^o |k|Ibkp)^. (4) 

Here Uk,bk,dk denote the Fourier modes of u,b and a 
respectively. The total dissipation of energy is given by 
Eg = e“ -I- e+, similarly for the square vector potential 
eA=e-+ e+ and vorticity = e" -h 

The injection rate for the enstrophy is dehned as R = 
{unpui) = luk'j and for the square vector potential as la = 
{a(j)a) = Rkj'^. Conservation laws then result to the 
relations = I a and = R + h- In the yf = 0 case 
we also have R = e . 

In the limit of large Reynolds number Re ^, Re~ —> oo 
all injected conserved quantities are transported in scale 
space from the forcing scale to the dissipation scales by a 
hux caused by the nonlinearity. Their dissipation rate is 
equal to the sum of their forward and inverse hux rates. 
The hux of kinetic energy H,^, enstrophy Hj,, total energy 
Hg and square vector potential H^ in the Fourier space 
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are given by, 

=(u<(u- Vu)) 

ng(fc) =(u^(u • Vu — b • Vb) + b^(u • Vb — b • Vu)) 
n^(fc) =(afc(u-Va)) 

nf,(fc) =(w^(u-Vw)) (5) 

where the notation represents the Fourier filtered field 
g so that only the modes satisfying the condition |k| < fc 
is being kept, see [T5] . 

Finally the distribution of energy among scales will 
be quantified through the kinetic and magnetic 
energy spectra: 


Euik) = 

fc<|q|<fc-|-l 

(6) 

Eb{k) = 

/c<|q|<fe-|-l 

(7) 


C. Numerical runs 

The aim of the present work is to reveal the mecha¬ 
nisms under which the cascade of the ideal invariants is 
changing as the control parameter /r/ is varied in the limit 
of large Reynolds numbers and box sizes. Due to compu¬ 
tational limitations we can not have all three Re'^ , Re~ 
and kfL large for the same runs. As an alternative two 
sets of runs have been performed. In the first for a mod¬ 
erate value of Re^ the parameter /i/ is varied for four 
different values of kfL and Re~ . The parameters for 
these sets of runs are shown in table HI Part of the data 
from these runs were presented in [5]. Here we present a 
more thorough analysis focusing on the mechanism that 
takes place during the observed transition. For complete¬ 
ness some of the findings in [5] are also shortly presented 
and discussed here. In the second set of runs /i/ was 
varied for a fixed value of kjL = 8 and Re~ and five 
different values of Re~^ to check the large Re'^ behavior. 
The parameters for these sets of runs are found in table 
[n| By comparing the different values of Re^ tested we 
can verify if the examined runs have reached a large Re^ 
asymptotic behavior. Finally, a third set of runs was per¬ 
formed with kfL = 4, i?e+ = 1.4 10 ^ and Re~ = 1.8 10 ^, 
with focus on the development of energy spectra in the 
small scales. 


III. CRITICAL BEHAVIOR 
A. Dimensional relations 

The relations between the injection rates, the forcing 
amplitudes and the amplitude of the fluctuations can 
be constrained by simple dimensional and scaling argu¬ 
ments, that we are using through out this paper. For 
large values of Re~^ the injection rates are expected to 


TABLE I: Numerical parameters of the first set of DNS 
keeping i?e+ = 1400 fixed. T expresses typical duration 
of the runs in units of 1 / y/fokf. 


Case 

A1 

A2 

A3 

A4 

kf 

8 

16 

32 

64 

N 

512 

1024 

2048 

4096 

Re~ 

2.3 10"^ 

7.410® 

2.310® 

7.6 lO’’ 

Re+ 

1400 

1400 

1400 

1400 

T 

2000 

600 

342 

300 


TABLE II: Numerical parameters of the second set of 
DNS keeping Re~ = 23000, kf L = 8 fixed. T expresses 
typical duration of the runs in units of 1 / f^kf. 


Case 

B1 

B2 

B3 

B4 

B5 

kf 

8 

8 

8 

8 

8 

N 

256 

512 

1024 

2048 

2048 

Re~ 

23000 

23000 

23000 

23000 

23000 

Re+ 

1.4 10^ 

1.410® 

1.410^ 

1.410® 

1.410® 

T 

5000 

2000 

600 

342 

300 


become independent of the dissipation coefficients in the 
small scales, and to depend only on the amplitude of the 
fluctuations at the forcing scale. We obtain then 

lu oc CpfoUf h oc Cpfiffobf. (8) 

Here by it/ and bf we denote the amplitude of the veloc¬ 
ity and magnetic fluctuations respectably at the forcing 
scale. In the presence of an inverse cascade the coeffi¬ 
cient Cp can depend on the amplitude of the large scale 
fluctuations that sweep the small scales and cause a fast 
decorrelation of the forcing and the velocity and mag¬ 
netic fluctuations that reduces the injection rate m- 
Thus in general is a function of Re but indepen¬ 
dent of i?e^. The injected energy, enstrophy and square 
vector potential is transferred to large or small scales by 
the nonlinearities. Assuming locality of interactions the 
fluxes and dissipation rates at small and large scales can 
be estimated to be 

= H^ oc u^/i 

= n| oc uebj£ (9) 

4 =n! 

where the relation for LI only holds in the /i/ 1 limit, 

and the relation should be replaced by ueiuj+cbyjt 
when magnetic effects become important. The propor¬ 
tionality coefficient in these relations is in general a non¬ 
trivial function oi fif and they can become zero in the 
case of critical transitions. This dependence on fif as 
well as the validity of the locality assumption can not be 
concluded from dimensional arguments and needs to be 
extracted from numerical simulations and this is one of 
the main objectives of this work. 

The dependence on /r/ however can be guessed in some 
limits. For example in the case that /r/ <C 1 the velocity 
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field is not affected by the magnetic field and thus the 
usual relations of 2D HD turbulence hold with ui oc H^^^ 
at small scales and being independent of /ij. In the 
same limit using the two relations and at the forcing 
scale we obtain, 


— fx flf, 


and ^ 

J-li. 


n 


ju4 

^Kf 






( 10 ) 


The balance of the fluxes with the small scale dissipa¬ 
tion rates also leads to the prediction of the dissipation 
length scales ii, that are given by £ 1 , oc £/ 
in the case the forward total energy cascade is dominant, 
oc if in the case the forward square vector 

potential cascade or the enstrophy cascade is dominant. 
Implications of these relations have been tested in the 
sections that follow. 


B. The limit of Re , kfL — > 00 


these two critical values no inverse cascade exists and 
both conserved quantities cascade only forward [8]. 



FIG. 2: e~/eE,e^/e^ are shown as functions of fj,^ for 
different values of Re~ as mentioned in table |l] Blue, 
continuous line corresponds to Ce/Ce and green, dashed 
line to e^/e^ with darker shades denote larger 
Re~, kf L as, mentioned in the legend. 




FIG. I: The normalized dissipation rates 
are shown as functions of p/ for different values of Re~ 
as mentioned in table |l] Blue, continuous line 
corresponds to (-~lf-E green, dashed line to e“/e^ 
with darker shades denote larger Re~, kf L as 
mentioned in the legend. 

We begin by investigating the large Re~ and large 
box size limit. Figure shows the normalized dissipa¬ 
tion rates e~feE, and as functions of /i/ for the 

four different values of kfL,Re~ from the set of runs 
described in table |T] Different symbols denote different 
values of kfL = (8,16,32,64) as mentioned in the leg¬ 
end. As the domain size and Re~ is increased the large 
scale energy dissipation c^/Ce goes to zero at a critical 
value oi fif implying that in the large box limit the tran¬ 
sition is critical: there is a value of the control parameter 
fif for which the amplitude of the inverse cascade be¬ 
comes exactly zero. Similarly the inverse cascade of the 
square vector potential goes to zero e“/e^ —t 0 in the 
large box limit at a different critical value of /i/. The 
two critical values of that are observed are denoted as 
Pec — 0.22 for which the inverse cascade of energy stops, 
and the other for which the inverse cascade of square 
vector potential starts denoted as p^c — 0.26 . Between 


FIG. 3: pe is shown as a function o^/ for different 
values of i?e“, kf L for runs in table [I Darker shades 
denote larger Re~, kf L as mentioned in the legend. 


As discussed in the previous section the alternative 
control parameter p^ can be used that is closer to theoret¬ 
ical models. The normalized dissipation rates as a func¬ 
tion of Pe are shown in hgure The shape of the tran¬ 
sition curve appears to change when the same data are 
presented as a function of p^. First the two critical points 
appear to merge to one, denoted now as pec — 0.045 . 
Second the two exponents appear that can be fitted close 
to the critical value are different from the p/ case and 
are both closer to one. These points imply a non-trival 
relation between p^ and p/. This relation is plotted in 
figure where the observable p^ is shown as a function 
of p/ for different values of Re~. Glose to the critical 
point fj-e becomes a non-monotonic function of p/. This 
non-monotonic behavior becomes more pronounced as we 
increase the domain size and is localized around the crit¬ 
ical points in the system. As a result for the critical 
value of pe we can find three values of p/ that satisfy 
Pe(p/) = Pec- Two of these values are the critical values 
Pec and p^c found in figure]^ and the third to a value 
in between. 
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The parameter /Xe is thus a parameter worth consid¬ 
ering for future studies for which the energy injection 
rate is fixed rather the forcing amplitude. However at 
the present study we will stay with description, since 
the uncertainty that comes from the measurement of /ie 
makes it difficult to interpret the results close to the crit¬ 
ical point. 


C. The limit of Re'^ — ^ oo 


The other limit considered is the limit of large Re~^ for 
fixed kf L = 8 and Re~ = 23000 (see table|^. A series of 
runs for five different values of Re'^ are made. For these 
set of runs the box size and Re~ are not large enough for a 
clear manifestation of criticality (that appears only in the 
large kf L and Re~ limit). For this reason we cannot dis¬ 
tinguish between the two critical points /i^c and that 
appear as one and we will simply refer to it as /ic. Figure 
4a shows the normalized dissipation rates at large scales 
as a function of /r/ for different Re'^. The critical point 
moves to smaller values of the control parameter as we in¬ 
crease Re'^. However, rescaling the control parameter f 


by multiplying with (Re 


+ sl/2n 


makes the curves collapse 


on top of each other as shown in Figure This is not 
so surprising since the magnetic field is being stretched 
and amplified at the smallest scales and its amplitude de¬ 
pends on the viscous cut-off length-scale. The rescaling 
factor can be derived assuming that the transition takes 
place when the magnetic field and the velocity field are 
of the same strength. From standard phenomenological 
arguments for constant flux (see section HI A) we have. 


2 ^ 1 
Of oc H^ — oc H^ 

ue nr 


( 11 ) 


The magnetic field at the viscous/ohmic scale £i, is then 
-1 /n \i/2 

given by, b^, oc a,y£^, o: Uf ( ) F^- Finally the 

small scale dissipation length scale £i, can be estimated 
by balancing the flux of A with its ohmic dissipation rate 
that leads to £,y oc £f Assuming that the 

transition takes place when bi, ^ Uf, and using equation 
[To] we obtain 


/Tc oc ^ oc (i?e+) (12) 

This leads to the observed scaling for the critical point 
/ic as a function of the forward Reynolds number. It 
is important here to note that in this scaling argument 
the transition occurs when there is an equality between 
the magnetic held at small scales and the velocity held 
at large scales. The transition is thus controlled by a 
nonlocal interaction between the small magnetic scales 
and the large velocity scales. 



(a) 



FIG. 4: The top plot shows e^/e^ and as 

functions of /ij for different values of Re'^ mentioned in 
the legend for runs in table jH] Dotted line represents 
e“/eg while the continuous represents e“/e^. Darker 
shades denote larger values of Re'^ as mentioned in the 
legend. The plot below shows the same data with a 
rescaled x axis /i/ (i?e+)^^^" 


D. Injection rates and fluctuations 

The relation between the injection rates and the am¬ 
plitude of fluctuations is the most basic relation in steady 
state turbulence. The behavior of these quantities as the 
control parameter approaches close to the critical point 
is of primary interest as they allow to estimate the rate 
energy cascades to the dissipation scales. Shown in flg- 
ure are the injection of kinetic and magnetic energies 
lu = = lak'f as functions of /r/ for runs 

in table ITj We note that fij was varied by keeping Fu 
hxed and changing the value of Ff,. Increasing values of 
Re~,kf L are denoted with increasing shade of color of 
the lines. Four values of Re~ are shown. The two verti¬ 
cal black lines at /r/ = 0.22, 0.26 denote the proximity to 
the critical points. 

In the limit /rj <C 1 the quantity reaches an asymp¬ 
totic value since the magnetic held plays a passive role 
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FIG. 5: Figure shows on top /„ and on bottom If,, the 
injection energies in the kinetic and the magnetic field 
as functions oi fj,f. Darker shades denote larger values 
of kf L, Re~ as mentioned in the legend for runs in 
tablel^ The black lines denote values = 0.22,0.26 
denoting proximity to the critical points ^c- 


FIG. 6: Figure shows on top uj and the bottom 6^, 
representing the square of the velocity field at the 
forcing length scale kf and twice the magnetic energy as 
functions oi fif. Darker shades denote larger values of 
kf L, Re~ as mentioned in the legend for runs in table 
The black lines denote values /i/ = 0.22, 0.26 denoting 
proximity to the critical points fj-c- 


and the flow is close to 2D HD. The value at the forc¬ 
ing scale is suppressed as kf L increases due to sweeping 
effect of the large scales from the inverse cascade [l9]. 
As expected the magnetic injection rate /{, a pproac hes 
zero in the /ry —^ 0 limit as Ib oc see section 


III A 


As 

we increase p,/ both injection rates /«,/{, increase. Glose 
to the first critical point fi^c both injection rates vary 
smoothly. Between the two critical points the injection 
rates sharply increase and at the second critical point 
/i^c both injection rates vary continuously but their first 
derivative with respect to fj,f appears to have a discontin¬ 
uous change. At large values of /ry the magnetic energy 
injection rate further increases while the kinetic energy 
injection rate is suppressed. Note however that for values 
of /iy shown in figure]^ Ib is always smaller than It 
is also worth pointing out that if instead of fj,f the data 
were plotted using both injection rates would appear 
discontinuous at the critical point. 


The injected energy at the forcing scale is transported 
in the dissipation scales by nonlinear interactions. As dis¬ 
cussed before the transition occurs when becomes of 
the same order as uj . Here is twice the total magnetic 
energy, 6^ = (|bp). Thus close to the transition the dom¬ 
inant role for the energy transport at the forcing scale is 
played by the velocity fluctuations at the same scale as 
the forcing and by the magnetic field fluctuations at the 
small viscous scales where the magnetic field is strong. 


The magnetic field at the forcing scale b^ is weaker than 
by a factor as discussed in section 


me 


^^ In 

figure® we plot the kinetic energy at the forcing scale 
(define as the energy of the Fourier modes with fc = |k| 
in the range [kf — l,kf + 1]) and the magnetic energy b^ 
that is dominant in the small scales. Uj becomes maxi¬ 
mum between the critical points where both the inverse 
cascade of energy and square vector potential vanishes, 
suggesting a pile up at the forcing scale kf. The magnetic 
energy b^ is much smaller than before the transition. 


b^ increases smoothly and rapidly after the first critical 
point. Glose to the second critical point there is a sharp 
change for both quantities, b^ shows a discontinuous first 
derivative while the derivative of uj: could also possibly 
diverge (from the right) as /ry is decreased to the 
value. 


E. Balance relations 

Having established the general behavior of the injec¬ 
tion rates and the amplitude of fluctuations close to the 
critical points we can attempt to give a phenomenological 
description of the cascade process. In 3D turbulence and 
in the limit of small dissipation coefficients one can argue 
by dimensional reasoning alone that the energy dissipa- 
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tion rate will be proportional to the cubic power of the 
fluctuation amplitude and inversely proportional to the 
forcing length scale Cg The non-dimensional 

proportionality coefficient can in principle depend on 
the forcing functional form. The same scaling holds for 
2D turbulence but for the dissipation rate in the large 
scales e“. 

The situation is different for the 2D MHD case, because 
there are two processes involved: vortex-shearing that 
amplifies energy in the large scales [20], and magnetic 
field line stretching that removes energy from the large 
scale vortices and amplifies energy in the small scales. 
Thus in principle we need to define two terms. The first 
one responsible for the inverse cascade is proportional to 
the cubic power of the fluctuations at the forcing scale 

kf. The second one that couples small scale magnetic 
field to the forcing scale shear is proportional to kf 
and suppresses the energy transport to the large scales. 
We thus expect for the inverse cascade: 

^~E=C~{u)kf-aufb^kf) (13) 

where C~ and a are non-dimensional constants that need 
to be determined from the data. As a proof of concept we 
plot in figure the quantities v?jkf/t~ and Ufb^kf /e~. 
For clarity only the case kf L = 16 is shown. Indeed these 
two curves show that close to the critical points the term 
responsible for the forward cascade Ufb^kf becomes of 
the same order as the term for the inverse cascade u^fkf. 
For values /r/ <C /r^c the term u^fkf dominates while for 
^ fijic the term Ufb'^kf is bigger. In figure we 



FIG. 7: The quantities u^fkf/e~ and Ufh'^kf/e~ are 
shown as functions of for the case kj L = 16, see 
table 1^ 

plot the ratio {u^fkf — aufb'^kf)/{u^fkf) for the different 
values of /i/ and Re~ of the runs presented in table 
The coefficient a has been chosen so that the numerator 
of this ratio becomes zero at the critical point. The data 
collapse to a single curve indicating that the functional 
form of this curve is independent of the large scale param¬ 
eters Re~ and kfL. We note however that the power-law 
behavior close to the critical point in figure [^differs from 
the one observed in figure This indicates that the co¬ 
efficient also depends on fif or alternatively mean 



FIG. 8: The ratio {u^fkf — aufb^kf)/{u^fkf) is shown as 
a function of /r/. Darker shades denote larger values of 
kfL, Re~ as mentioned in the legend for runs in table 

square values (such as b^,uj) cannot provide an accurate 
description near the transition and higher order statistics 
need to be taken into account. 

The case for the inverse cascade of square vector po¬ 
tential is different because it is not as easy to identify 
the different processes involved. The simplest possible 
dimensionally acceptable form for the dependence of the 
inverse cascade on the amplitude of the fluctuations is 
given by: 

e; = C-Ufa)kf = C-Ufb}kJ^ (14) 

where now C~ is a function of the control parameter 
^f. Plotting the ratio e~/ufb'jkj^ for all the data of the 
runs in table|||indeed collapses the curves to a single one, 
shown in figure This indicates that the transition to 
the inverse cascade of the square vector potential can be 
expressed by a “Kolmogorov-coefficient” that varies with 
fj,f and becomes zero at the critical point. Note however 
that this simplification originates from our lack of knowl¬ 
edge of the precise mechanisms involved in the transition 
for square vector potential. In principle even this tran¬ 
sition could result from the competition of two processes 
(one forward e.g. by eddy stretching and one inverse e.g. 
by a negative turbulent diffusivity) that equate at a crit¬ 
ical value oi fj,f. However, at present we cannot identify 
these processes from the attained data, and leave such 
considerations for future work. 


IV. IDEAL INVARIANTS ACROSS SCALES 
A. Flux of conserved quantities 

The flux of a conserved quantity across scales gives 
the rate that the nonlinearity transports this conserved 
quantity across scale space. Figure shows the four 
fluxes , Hj^, Hg, as a function of the wave-number 
k for different control parameter fif for the kfL = 8 case, 
corresponds to case B3 in tablejnj We start with the case 
Hf = 0.0001 <C 1, where the vector potential acts like a 
passive scalar. An inverse flux of kinetic energy H^^ < 0 is 








FIG. 9: The Kolmogorov constant C~ is shown as a 
function of the control parameter for different values 
of Re ~, see table |lj Vertical black lines are 
= 0.22, 0.26 denoting proximity to the critical points 
fif. Darker shades denote larger values of kf L, Re~ as 
mentioned in the legend. 


observed at wave-numbers smaller than kf and a forward 
flux of enstrophy > 0 ioi k > kf. For k > kf the flux 
of square vector potential > 0 has a forward cascade 
as expected for passive scalars. Since the magnetic field 
is small the flux of total energy is equal to the flux of the 
kinetic energy with a correction of the order fcj, = 
+0{fL‘f). For larger values of the conservation of ft 
breaks down progressively and the vector potential stops 
acting like a passive scalar. The Lorentz force injects 
enstrophy into the flow leading to non constant flux of 
IIj^, The enstrophy flux is modified most strongly at 
the largest wave-numbers where the curl of the Lorentz 
force Fl = e^j x b is larger. For values close to the 
critical point ^Ji,f we have a simultaneous forward 

and inverse cascade of total energy. As /r/ becomes larger 
than fL(. the cascade of becomes strictly forward and 
an inverse cascade of is established. 

To further demonstrate the two processes involved in 
transferring energy to and from the large scales we calcu¬ 
late = rig — Ilg. Ilg expresses the energy flux from 
terms that do not involve the magnetic field, while 11^ 
is the energy flux from terms that involves the magnetic 
field. Figure [TT] shows 11^, Ilg, Ilg as a function of k for a 
fjLf which is close but smaller than the critical point. The 
two components of the fluxes 11^ and 11^ are not constant 
neither in the forward nor in the inverse cascade. Ilg is 
negative (inverse) at almost all wave-numbers while 11^ 
is positive (forward) at all wave-numbers. Thus energy 
is transported by the velocity field from the small scales 
to the large scales, while the opposite role is played by 
the magnetic field that transports energy from the large 
scales to the small scales. The weak but constant en¬ 
ergy flux observed at large scales is then the result of a 
counter balance of the inverse flux of the kinetic energy 
with the forward cascade of magnetic energy. This in 
part also explains the large fluctuations observed in the 
flux of energy in [S]. 



FIG. 10: The four plots from top are normalized fluxes 
of total energy Ilg/eg, square vector potential 11^/e^, 
kinetic energy 11^ /, and the enstrophy IIj, /as 
functions of wave-number k for different values of the 
control parameter fif. Darker shades denote larger 
values oi fif as mentioned in the legend. 


B. Spectras 

In the previous work [8], the spectra of kinetic and 
magnetic energy £'„(fc), Eb{k) were shown near the tran¬ 
sition focusing on the behavior of the large scales. For 
values oi fif smaller than the critical value the in- 
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FIG. 11: The normalized flux of 11^,, II^, 11^ for the 
control parameter /r/ < /ic is shown as a function of k, 
with parameters corresponding to case B3 in table 


verse flux of kinetic energy formed spectra proportional 
to at large scales for the velocity field. The mag¬ 

netic field showed equipartition of A among large scales 
which led to the spectra Eb(k) ~ For values of /r/ 
larger than the magnetic energy starts to form an ex¬ 
ponent of k~^^^ due to the inverse cascade of A mm- 
Here we focus on the spectra at wave-numbers larger 
than kf and how they change as the parameter /i/ is 
increased from infinitesimal values. In order to maximize 
the inertial range for the forward cascade we performed 
runs with kf = 4, with Re'^ = 1.410^,i?e“ = 1.810^. 
The kinetic energy Eu{k) and the magnetic energy Eb{k) 
spectra are shown in figure[^for few values of the control 
parameter /i f . For these set of runs the parameter /i f was 
kept smaller than its critical value /i^c as the spectra 
in the forward cascade did not change significantly for 

^ A^ec- 



FIG. 13: The ratio of the dissipation at small scales of 
the magnetic and the kinetic field /e^ are shown as 

function of the rescaled control parameter fif 
Four different value of Re'^ are considered as mentioned 
in the legend, for runs in table ITTj The continuous black 
line represents the value /ej = 1- The dotted black 
line shows the scaling at the passive scalar limit. Darker 
shades correspond to larger i?e+. 


For the lowest value of /r/ the magnetic energy spec¬ 
trum is smaller than the kinetic energy spectrum at all 
scales and corresponds to the passive scalar limit. The 
forward cascade of enstrophy leads to a kinetic energy 
spectrum Eu oc The magnetic field on the 

other hand is proportional to the gradients of a passively 
advected scalar, forms a spectrum Eb oc 
However since the magnetic spectrum has a larger expo¬ 
nent than the kinetic energy spectrum the two spectra 
will meet at a wave-number k^ that can be obtained by 
equating the two relations for E^ and Eb and leads to 



FIG. 12: The spectra of the kinetic E^ and the 
magnetic energy Eb as functions of k are shown for 
different value of ftf. Dotted line represents the 
magnetic energy spectra and the continuous line 
represents the kinetic energy spectra with darker shades 
corresponding to larger fif. Black lines represent the 
possible exponents of the spectra. 


kra (X (X fj,/^'^kf. (15) 

This however is realized provided that the inertial range, 
limited by the dissipation, is long enough so that k^ < 
ki,. This in turn implies that the transition from the 
passive scalar regime {b^, <C u^) to a regime at which 
magnetic and kinetic energy are comparable at the small 
scales (6i/ ^ u^) occurs for a value /i/ that we will refer 
to as To have an estimate for we look at the 
forward dissipation length scale in the passive limit given 
by in = k~^ oc kj^ By equating kj^ with k^, 

we have the scaling, 

k-NL ^ ~ (16) 

This scaling is clearly demonstrated in figure where 
the ratio /e+ is plotted as a function of fif {Re'^y^^ 
that collapses all curves together for different values of 
Re^ from Table ^ The ratio e)J“/e+ expresses the ra¬ 
tio of magnetic tolcinetic energy weighted at the small 
scales. Glearly for all runs /e^ is smaller than unity 
for all values of /r/ such that iJ,f{Re'^y^^ < 3 that marks 
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the beginning of the nonlinear feed back of the magnetic 
field. Note that the onset of the nonlinear behavior 
is different from the the critical value /Xc oc 
that marks the end of the inverse cascade. 

For values of larger than the passively ad- 
vected magnetic field spectra appears to still hold but 
only for wave-numbers for which k < km- This can be 
seen in figure for the kinetic and magnetic spectra 
in the range of wave numbers for which Ef, <?; For 
larger wave numbers the observed spectral slopes change. 
For sufficiently large /i/ and for values k > km there 
is a new power law behavior observed that is close to 
^-5/3 fQj. both Eh and E^ spectra. This is in agreement 
with the 2D MHD prediction of a constant forward en¬ 
ergy flux to small scales. Thus the initial power laws 
Eh{k < km) ^ k~^^,Eu{k < km) ~ transitions to 
Eh{k > km) k~^/^,Eu{k > km) k~^/^ at a fixed 
value oi fif. As we increase the value of the magnetic 
energy becomes larger thus the point km moves closer 
to the forcing length scale kf as expected from the es¬ 
timate in equation |15[ These arguments finally break 
down when /r/ comes close to the critical value /i^c- 


V. SPATIAL BEHAVIOR 


A. Structures 


The different processes and the different phases dis¬ 
cussed in the previous sections have a direct impact on 
the structures formed. For values of /xy < the flow 
forms large scale vortexes due to the inverse cascade of 
energy and weak small scale vortex filaments due to the 
forward cascade of enstrophy. Properties of these struc¬ 
tures have been extensively discussed in the literature 
|20| . At the same time the vector potential, that acts as a 
passive scalar, forms also filamentary structures p51124j . 
This behavior starts to change as fif becomes larger than 
/Xjy^. Figure 14 shows the vorticity w and the current j 
for /x„j^ < /^/ < and kfL = 4. The large scales struc¬ 
tures seen in w is due to the strong inverse cascade of the 
kinetic energy similar to 2D HD systems. The small scale 
structures in uj field however are not similar to the ones 
of 2D HD turbulence but closer to 2D MHD systems. 
The current density in the same figure leads to the same 
conclusion. The current is concentrated in thin filaments 
aligned with the large scale shear as for the case of an 
advected passive scalar. However these filaments are not 
‘straight elongated structures’ as the case of the gradi¬ 
ents of a passive scalar but rather the filaments ‘wiggle’ 
varying along the direction of the shear, probably due to 
the appearance of nonlinear MHD (Alfven) waves. 

Figure [TS] shows the kinetic and the magnetic energy 
Eu,Eh for a value fif < /x^c and kfL = 64. A large 
portion of kinetic energy is in small scales with the ex¬ 
istence however of scales much larger than the forcing 
scale, a clear evidence of an inverse cascade. The large 
scale structures in the kinetic energy are not space filling 



FIG. 14: Contour plots of the vorticity uj on top and 
the current j in bottom are shown for a value 
Mjvi < M/ < Me with kf L = A. Both large scale and 
small scales structures are found due to the dual 
cascade of energy and square vector potential. 


suggesting that the inverse cascade is weaker than the 
pure 2D HD. Magnetic energy is concentrated only in 
the small scales. However we note that there is an anti¬ 
correlation between the kinetic field and the magnetic 
field: Regions of large velocity structures are correlated 
with regions of weak magnetic field and vice versa. Thus 
it appears that some regions of the flow act as 2D HD 
cascading energy inversely and other regions act as 2D 
MHD cascading energy to the small scales in accordance 
with the process suggested in figure El 

Finally figure shows the magnetic energy Eh and 
the square vector potential for a value fif ^ fi^c and 
kfL = 64. The kinetic energy (not shown) does not 
show any large scale structure. The magnetic energy is 
clearly small scale with a filamentary structure. However 
the current filaments are arranged so that the vector po¬ 
tential forms large scale islands. 
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FIG. 15: Contour plots of the kinetic energy field on 
top and the magnetic energy field in bottom are shown 
for a value /r/ < with kf L = 64. 


B. Spatial statistics 

To be able to quantify the qualitative description 
given in the previous section we calculate distribution 
of differences of velocity 6uj^ (x, r, t) and vector poten¬ 
tial 6a (x, r, t) between two points x and x -|- r. Here 
subscript L denotes longitudinal component of the vec¬ 
tor 6u^ = i5u^ • f . Studies of 2D HD [Ml [251 ES] have 
suggested that the inverse cascade in 2D turbulence is 
self-similar (and possibly conformal invariant [27j l in the 
sense that the pdfs of velocity difference from two points 
at distance r (that lies in the inverse inertial range) take 
the form oc f(5u^/r°‘) for some exponent a . 

The probability distribution function (pdf) of 
the normalized longitudinal velocity difference 

(x, r, t) / (x, r, t)) and the normalized 

1 1 ^ 

vector potential difference da (x, r, t) / (x, r, t)) 

for different values of r are shown in figure [TtI The pdfs 
were calculated using the results from the resolution 



FIG. 16: Gontour plots of the magnetic energy field on 
top and the vector potential in bottom are shown for a 
value /T/ ^ with kf L = 64. 

runs in table |l] Different values of the control parameter 
fjLj are shown by successive vertical shifts in the y-axis. 
Different colors corresponding to different values of r 
ranging from r~^ = 0.9 x kf to 0.1 x kf corresponds to 
the inverse cascade range. Black curves correspond to a 
normalized Gaussian curve of unit variance. For all the 
values of r the pdfs collapse on each other indicating 
self-similarity. This is seen to be consistent across 
transition as the value of /i/ varies from 0.05 — 0.5 with 
He ~ 0.23 and kf L = 16 see table [Tj Similar behavior is 
observed for larger values of kf L = 32,64. For pdf of 
close to the transition the presence of larger tails 
suggests higher probability of extreme events. 

From the different moments that we can examine, 
of particular interest are the third order moments that 
are related for homogeneous and isotropic flows by 
Karman-Howarth type relations for the total energy. The 
Karman-Howarth relation for E is given by [55], 

- = (Sul (x, r, t)) - 6 {bl (x, t) 6u^ (x, r, t)}l7) 
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For the energy balance we have a competition be¬ 
tween the inverse cascade of kinetic energy created 
by (x, r,<)) and the forward cascade created by 
(x, r,t)6^ (x,r,<)). Thus the direction of cascade 
is determined by which term dominates. 




FIG. 17: The normalized pdfs - 
5u^ (x, r,t) / (x, and 

5a (x, r, t) / (^5a? (x, r, t)) on top and bottom 
respectively are shown for different values of 
fif = 0.05,0.21,0.26,0.5 represented on the side. 
Different values of are shown by displacing the pdfs 
vertically. Different colors represent different values of 
r. Black lines represent the gaussian curve with unity 
variance. 

The pdfs of 6 (6^ (x, r, t) (x, r, t)) and 
(x, r, t)) slightly before and after transition 
are shown in figure The values of the control 

parameter are /.t/ = 0.21 and 0.26. Self-similarity across 
all r in the inverse cascade scales allows us to average 
across r in the inverse cascade inertial range. The 
symmetric part of the pdf (where for a function f{x) it is 
defined as (/(a;) -I- f{—x)) /2) is shown by the black line. 
The asymmetry in the curve gives an indication for the 
direction of the cascade of the energy across scales. For 
(top panel) the pdfs have exponential asym¬ 
metric tails with the asymmetry in favour of the positive 



FIG. 18: Pdfs 5u^ in red and 6b^6u, in blue are shown 

L L 

for two different values of fj,f. fj,f = 0.21 < fic is the one 
on top and /i/ = 0.26 > /ic is the bottom one. Black 
lines denote the symmetric part of the pdfs. 


values. Thus based on the Karman-Howarth relation the 
(x, r,t)) term leads to a flux of energy towards the 
large scales and the term —6 (6^ (x, r, t) (x, r, t)) 
leads to a flux of energy towards the small scales. For 
values of larger than the critical point p./ > the 
tails remain non-Gaussian but the asymmetry in the 
pdfs of both terms is diminished, thus the exchange of 
energy with the large scales is suppressed. 


VI. SUMMARY AND PERSPECTIVES 


We have studied by an extensive number of numer¬ 
ical simulations the transition from 2D HD turbulence 
to 2D MHD turbulence, varying systematically the mag¬ 
netic forcing and the involved Reynolds numbers that 
allowed us to identify different phases of the turbulent 
state of 2D MHD. Our findings are summarized in fig- 
The top panel shows the variation of e^/e^ 
“/e^ as a function of the control parameter /i/ 
for fixed Re^ >> l,kf L ^ 1. For small values of 
the first phase of 2D MHD is met where the mag¬ 
netic field is advected passively without altering the hy- 


ure 

and 
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FIG. 19: Top panel: the variation f -~!^ 
function of Denoted are the critical points ^IectI^ac 
and the region where the nonlinear behaviour starts 
denoting the Reynolds number calculated 
based on the order of laplacian n. Lower panel: The 
phase space diagram in the parameter space Re'^, fif. 
Re~, kf L is assumed to be large or infinite. 


drodynamics. Thus energy cascades inversely and the 
square vector potential cascades forward. The magnetic 
field becomes active first in the small scales when /r/ 

is larger than ex (Re"*") In this phase, part 

of the energy cascades forward to the small scales and 
part of the energy still cascades inversely to the large 
scales, while the cascade of the square vector potential 
remains forward. Further increasing we reach the 

critical point oc {Re'^) where the inverse en¬ 

ergy cascade stops and all energy cascades forward. For 


slightly larger values we meet the second critical point 
that also scales as {Re~^) and marks the be¬ 
ginning of the inverse cascade of A. Between these two 
points /ijjc < M/ < k'AC both fluxes are strictly forward. 
For values of /i/ larger than the phase of 2D MHD 
turbulence exists where the energy cascade is forward and 
part of the square vector potential cascades inversely and 
part forward. As tends to very large values we ex¬ 
pect that forward cascade of A will tend asymptotically 
to zero and A cascades strictly inversely. The different 
phases are shown in the Uf and Re~^ parameter space in 
the lower panel of figure [T^ In the phase diagram Re~ 
and kfL are assumed to be very large or infinity so that 
the transition clearly exhibits critical behavior with well 
defined critical values of /i/. The three different phases 
of turbulence are marked on this diagram and are sepa¬ 
rated by power laws of the form Hf oc (i?e+)^. Whether 
/igc and converge to a single point as i?e+ tends to 
infinity remains to be seen in future investigations. 

Our work also indicates that close to the two critical 
points the dependence of the injection rates and the fluc¬ 
tuation amplitudes do not necessarily have a smooth be¬ 
havior allowing the possibility that their derivatives with 
respect to /r/ to diverge or be noncontinuous. This is 
not unusual of course in phase transitions where differ¬ 
ent moments of the fluctuations, or susceptibilities scale 
as non-integer power-laws with the deviation from crit¬ 
icality [53]. To test such a prospect and measure these 
exponents precisely, many runs close to the critical points 
with both high Re~^ and Re~ are needed. Unfortunately 
current computational limitations restrict us from per¬ 
forming a study in such detail. 

Finally our study indicated that criticality arose from 
the balance of two counter-cascading processes one driven 
by the hydrodynamics (inverse) and one driven by the 
magnetic field (forward). Such counter-cascading mech¬ 
anisms also exist in other systems exhibiting variable in¬ 
verse cascade such as rotating and stratified flows 110 
It remains to be seen if these mechanisms also lead to 
critical behavior. 
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